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Abstract 



^ • The notion of the integral over the anticonimuting Grassmann variables is applied to 

^ \ analyze the fermionic structure of the 2D Ising model with quenched site dilution. In 

the A''-replica scheme, the model is explicitly reformulated as a theory of interacting 
fermions on a lattice. For weak dilution, the continuum-limit approximation implies 
the log- log singularity in the specific heat near Tc- |j 



> 

o 

^ 1 Introduction 

\D 

O 

Qs^ ' The site dilution provides, probably, the simplest way to realize quenched disorder in real 

On ■ magnetic materials. In this case, some amount of the magnetic atoms in a sample, chosen 

at random, are replaced by nonmagnetic impurity atoms. Another sort of disorder is bond 
dilution, in which case some of the lattice bonds are assumed to be broken. The two- 
dimensional Ising model (2DIM) is a natural object to analyze the effects of disorder on 

r~| , ferromagnetic phase transition since the exact analytic solution is known for this model 

in the pure case |j|]. The disordered 2D Ising model already has been analyzed during 
the last decades both by theoretical tools and in the precise Monte-Carlo experiments 
[^, ^, ^, m, p. The theoretical studies were merely concerned with the case of bond 

Cd \ dilution, making use of the fermionic interpretation of the problem |P, ^, ^. In this 

report, we apply a new anticommuting path integral technique 0, ^ ^ to clarify the 



fermionic structure of the 2DIM with quenched site dilution [10|. At the first stage, the 
partition function with fixed dilution is transformed into a fermionic Gaussian integral. 
The averaging over the disorder within the A^-replica scheme then results a lattice fermionic 
theory with interaction. Unexpectedly, the interaction in the exact lattice theory appears 
to be of order 2A'^ in fermions, where A^ — > is the number of replicas [^. An effective 
four-fermion interaction arises, however, in a continuum-limit field theory responsible for 
the critical behaviour near Tc for weak dilution due to the interplay of the short-wave 
and long-wave lattice fermionic modes as we pass to the continuum limit. In particular. 



^ Published in "Path Integrals from peV to TeV: 50 Years after Feynman's Paper", Proceedings of 
the 6th International Conference on Path Integrals from peV to TeV, August 25-29, 1998, Florence, Italy. 
Edited by R. Casalb uoni, R. Giachot ti, V. Tognetti, R. Vaia and P. Verrucchi (World Scientific, Singapore, 



1999) p.p. 137-141. hep-th/9906107 



this implies the double-logarithmic singularity in the specific heat near T^ for weak site 
dilution. The question about the effects of strong and moderated dilution, however, still 
remains to be an important open problem in the 2D Ising model, this holds true for both 
site and bond dilution. 

2 The model 

Let us start with a fixed distribution of the nonmagnetic sites over a lattice (fixed site 
dilution). The Ising spins, amn = ili are located at lattice sites, mn, with m,n = 
1,2,...,L running in the horizontal and vertical directions, respectively. Here L is the 
lattice length, L — > oo. To introduce site dilution, we accompany each Ising spin by a 
random variable ymn = 0, 1. The hamiltonian is [p!o| ]: 

— PH {y\a} = 2_^[bl ymnym+lnO'mnO'm+ln + ^2 ymnUmn+lO'mnC^mn+l] , (1) 

mn 

where 61^2 = PJi,2 suce the bond coupling parameters, Ji^2 sue the exchange energies, 
P = 1/kT is the inverse temperature. We assume ferromagnetic case: 61,2 >0. Noting the 
identity for a typical bond weight: exp ibyy'aa') = cosh.{b yy') + yy'aa' smh{byy'), which 
readily follows from aa' = ±1 and yy' = 0, 1, the partition function can be written in the 
form: Z {y} = R{y} Q {y}, where R{y} is a nonsingular spin-independent prefactor, to 
be ignored in what follows, and Q{y} is the reduced partition function | p!o| : 

Q{y} = Sp II (1 + ii l/mnj/m+ln0"mn<7m+ln)(l + i2 2/mnymn+l<7mn'7mn+l) , (2) 

y'^) mn 

where ti^2 = tanh 61^2 and a properly normalized spin averaging is assumed. Since we 
are interesting in quenched disorder, we have to average the free energy rather than the 
partition function itself. The standard device to avoid the averaging of the logarithm is 
the replica trick: 

-/3/Q{y} = lnQ{y}= hm l(Q^{y}-l). (3) 

In this scheme, we take A^ identical copies of the original partition function and average 
Q^ {y} ■ The formal limit A^ — > to be performed at final stages. In what follows, we 
assume the simplest distribution of the impurities in the definition of the averaging: 

W (ymn) =p6{l- ymn) + {I - p) 5 (ymn) , (4) 

where 5{ ) are the correspondent Kronecker's symbols, p is the probability that any given 
site, chosen at random, is occupied by the normal Ising spin, while 1 — p is the probabilty 
that the given site is dilute. 

3 Fermionic interpretation 

The partition function with fixed disorder (|2|) can be transformed into a Gaussian fermionic 
integral g following the method of the mirror-ordered factorization for the density ma- 

^ The rules of the integration over the purely anticommuting (Grassmann) variables were first intro- 
duced by F.A. Berezin |7|. In a close analogy with the bosonic case, a Gaussian fermionic integral of any 
kind can be expressed in terms of the determinant of the associated matrix |7| . For a short comment about 
the Gaussian fermionic integrals also see [jnl, [L2[. 



trix.g g |. This results the integral |To|: 



QlUs — / I I "■^mn"-^mn"-^mn"''^rnn ^'^P / , Y^mn^mn ~r "mn'^mn ' \^) 

mn mn 



amn 



^mn "mn 



where amn-,o,^n->^mn-,bmn ^^^ Grassmann variables. In turn, integrating out a part of 
fermionic variables from (S), we obtain the reduced integral |1C]: 



V 12/ J / i i rnn"'^mn Vmn 6Xp y ^ I Vran ^mn^mn ~r 

mra mra 

' \Pmn ~r ^mn) V'l^m—ln ''S'^mra— iJ ymn'^l'^2 C-m—ln^rnn—l i [Pj 

where Cmra,Cm,n are Grassmann variables, y^„ exp (y;;;^CmnCmn) = Umn + CmnCmn- The 
integrals (^) and @ are completely equivalent to each other and to (|2|). The advantage 
of the reduced representation like @ is that the corresponding action, in the pure case, 
explicitly illuminates the Majorana-Dirac structures of 2DIM already at the lattice level 







The averaging of (|6D over the disorder within the A^-replica scheme, see (jS]) and 
results the theory with interaction of the following kind pO[ : 



Q^{y} = / n n dci"r}dcj^2 n [p n e""™ + (i -p) n ^^^~^i^^ ] (7) 

mra a=l mn a=l a=l 

N N -, N 

P / 1 i 1 i nin "'^mn ^-^P / y I / ^ '-'mra ' ~ _l J_ ^mn ^mn ^ J 

mra Q!=l mra a=l " a=l 

where Smn is the replicated Gaussian action from (y) for the pure case, with y^n = 1 at 
all sites, and Amn = —Smn- Taking into account the nilpotent property of fermions in 
prefactor before the exponential, another possible choice is: Amn = ^1^2 Cm-ira^mn-i 10 ■ 
The continuum-limit field theory for weak dilution that follows from (^ is commented 
in the next section. A relatively simple form of interaction in (^ also provides grounds 
to try the cases of strong and moderated dilution, starting directly with the exact lattice 
integral (^). The analysis of this kind, however, have not yet been performed. 

4 The Gross-Neveu model (A^ -^ 0) 

To extract the effective continuum-limit field theory from (1^, we have to distinguish 
explicitly the higher- and low-momentum lattice fermionic modes (Fourier-harmonics) in 
the exact lattice action (^. Integrating out the higher-momentum modes in the first 

^ The starting point of the method is the fermionic factorization of the local spin-variable bond 
Boltzmann weights, {a\(j'). Introducing intermediate fermionic variables, we write (cr | a) {a\a'). After 
elimination of the spin variables, it is then possible to express the partition function in terms of the purely 
fermionic factors, (a j a'). In fact, the mixed factors like {a \a), arising by inserting the fermions, are neither 
commuting nor anticommuting with each other. So, one have to invent a special ordering for these factors, 
in their global products, in order the elimination of spin variables be really possible ISl 0] . In particular, 
the ordering problem is an obstacle to deal with the 3D Ising model. In two dimensions, however, the 
method is simple [H, M . Following this way, we come to ( 
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order of perturbation theory (weak dilution), we obtain an effective action for the low- 
momentum fields (large distances). The effective action appears, finally, in the form of 
the A^-colored Gross-Neveu model (A^ — > 0) ||lC 



5C~N = / d^ Y. [^N ^i"V^"^ + I V'i"^ (5l + Z d2) V!"^ (8) 

+ ^ 4"^ (-5i + ^ d,) 4"T + SAT [ E V'i^VfT } , 



a=l 



""^ v/2(iit2)c I' (^) v/2(tit2). 



P (^)' 4(tit2)c' 

where ipi,tp2 are the anticommuting Majorana components, ttTat and gfTv are the effective 
mass and charge, respectively. Here (A) and {B) are some lattice fermionic averages 
(short distances) explicitly calculated in [^]. The Gaussian part in (^) is the replicated 
Majorana action, corresponding to the pure case, with the mass term modified by disorder. 
The A^ = G-N model similar to (^) already has been analyzed intensively by DD-SSL 
as an effective theory near T^ for weak bond dilution |^, ^ ||, 13 1. The predictions are 



the double-logarithmic singularity in the specific heat and the logarithmic corrections in 
other thermodynamic functions. For more details see [Q, |5|, 0, |l^. The present analysis 
thus confirms the idea of the universality for small amount of impurities in the 2D Ising 
ferromagnets, for weak dilution. The situation is less evident, however, for strong and 
moderated site and/or bond dilution [^, |l^ |l^] . 

5 Added note 

Let / (ci,C2, ..., Cat) be any function of Grassmann variables ci,C2, ..,CAr. From the basic 
rules of fermionic integration, it follows: 

dcN---dc2dci f (Aici, A2C2, ..., Xncn) 
= \iX2... Xn dcN..-dc2dcif {ci,C2,...,cn) ■ (9) 

Infinite series of related identities also follows by differentiating (y) with respect to the 
parameters Ai, A2, ..., Aat. In particular, applying these ideas to the integral (0), we obtain 
the identity: 

r 2 / ^ \ '-^mn ^^ p -l-l ^mn'-mn ^ / ^ ) K-^'-'/ 

mn P=l 

where the averaging is assumed to be taken with respect to the non-Gaussian fermionic 
measure from (0). This identity can be used to check the consistency of approximations 
of any kind when dealing with lattice theory (1^). In a sense, the equation (jTH) can be 
viewed as the analog of ij? j^rn) = 9/2, 9 = kT, for the kinetic energy of a particle in 
classic Maxwell gas (one degree of freedom). In any case, the later identity readily follows 
by exactly the same method. 
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